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ABSTRACT 


An  experimental  and  theoretical  study  has  been 
undertaken  which  considers  large  deflections  of  a  light 
arch  under  the  action  of  a  radial  centrally  applied  load. 
The  arch  has  a  circular  span  and  the  ends  of  the  arch  are 
built  in.  Symmetrical  and  asymmetrical  buckling  is 
considered • 

Experiments  are  described  in  which  a  variety 
of  arch  dimensions  are  considered,  and  the  results  are 
compared  with  those  given  by  two  recent  theories.  These 
theories  are  also  compared  to  each  other. 

Discrepancies  arise  between  the  experimental 
considerations  and  the  theoretical  predictions.  It  is 
suggested  that  the  weight  of  the  arches  is  responsible 
for  such  discrepancies. 
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NOTATION 


Non-dimensional  amplitude  of  symmetrical  deflection 
where  a  =  A 


R£2 

Non-dimensional  deflection  corresponding  to  upper 
buckling-  load 

Non-dimensional  deflection  corresponding  to  lower 
buckling  load 

Non-dimensional  deflection  corresponding  to 
transitional  buckling  load 

Function,  where  a  =  d  -  4 

Width  of  arch  cross-section,  inches 

O 

Function,  where  d  =  n  a 
Strain 

o 

Function,  where  hn  =na 

1  ~T$u 


Function,  where  h2  =  3./  d  -  l\ 

u  \16  3  / 

Change  of  curvature  of  arch 
Number  of  terms  in  a  series 

Number  of  terms  in  a  series 

* 

Function,  where  p  =  4&uP  -  4 
Refers  to  the  rth  term  in  a  series 
Thickness  of  arch,  inches 
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NOTATION 


u 

v 

w 

X 

A 

B 

C 

Cm 

c2n 

c3n 

cAn 

% 

d3 

E 

H 

H 

I 


Function,  where  u  =  X 

7T2 

Tangential  displacement  at  middle  surface  of  arch 
Radial  displacement  at  middle  surface  of  arch 
Function,  where  x  =  oc 

J3 

Amplitude  of  deflection,  inches 
Amplitude  of  deflection,  inches 
Constant  of  integration 

"  I 

Definite  integral,  where  C-j_n  s  wn  dx 

Definite  integral,  where  C2n  = 


Definite  integral,  where  C3n  = 


Function,  where 
Summation,  where 

Summation,  where 
Summation,  where 


\n  =  ^wn)T  _ 


dx 

\dx  ) 

,  2  ,2 

dx 

\dx2  / 

x  =  0 

Di 


' C2n  C4n 

(C3n  -  KC2n)2 


D 


r 

T 

I 

D3  =  r 


cln  c3n  c4n 
(C3n  -  KC2n)2 

Cln  (2C3n  -  KC2n) 


I  L 


(C 


3n 


-  KC2n)2 


Young1 s  Modulus  (Appendix  A) 

Total  energy  of  system 

* 

Function,  where  H  =  HX 

Second  moment  of  area  of  arch  cross-section 
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NOTATION 


Constant,  where 

K  =  2.0246tt2 

Function,  where 

Kc  =  id 

36 4- 3u2  (9a  -  116) 

Tt 

it 

K1  - 

-12  +  3u2  (6  a2-  120)" 

tt 

tt 

K2  a  a 

-20  +  3u2  (-5T2  -  20) 

Applied  load,  lbs. 


Intermediate  buckling  load ,  or  energy  load 
Lower  buckling  load 
Transitional  buckling  load 
Upper  buckling  load 

Non-dimensional  load,  P  =  PR 

it2b]3 

o 

Non-dimensional  load ,  P  =  PR 

SI 

Radius  of  curvature  of  arch  centre  line 
Non-dimensional  strain  energy,  due  to  bending 
Non-dimensional  strain  energy,  due  to  axial  deformation. 
Non-dimensional  potential  energy 

2 

Non-dimensional  parameter,  where  X  =  JQj 

t 

Angular  coordinate  of  arch  (measured  from  centre  line) 
Half-angle  subtended  by  arch 

Non-dimensional  central  deflection,  where  c)=  A 

R 

Initial  state  prior  to  buckling 

Change  in  quantities  during  transition  of  arch  from 
unbuckled  to  buckled  state. 
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INTRODUCTION 

Little  work  considering  the  phenomenon  of  clamped 
arch  buckling  was  published  until  the  completion  of  the 
Second  World  War,  when  A.  van  Wijngaarden^  presented  a 
theoretical  treatise  upon  large  distortions  of  circular 
rings  under  the  action  of  a  concentrated  radial  load.  In 
his  paper,  by  introducing  boundary  conditions  to  adapt  to 
the  case  of  a  clamped  circular  arch,  he  derived  a  non-linear 
differential  equation  which  satisfied  the  geometrical 
conditions  and  the  equilibrium  state  of  the  arch.  This 
was  subsequently  solved  by  algebraic  transformation  and 
the  introduction  of  elliptic  integral  theory.  Symmetrical 
deflections  only  were  discussed. 

Experimental  investigation  was  also  undertaken, 
some  results  of  which  were  presented,  although  no 
description  of  this  work  was  offered. 

J.  S.  Kennedy2  reviewed  van  Wi jngaarden * s  paper 
much  later  in  1955  and  conducted  a  series  of  experiments 
with  light  thin  arches  whereby  he  induced  the  arches  to 
buckle  symmetrically  by  clamping  the  centre  of  the  arch. 

The  arch  centre,  therefore,  was  allowed  to  move  only  in  a 
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. 
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vertical  direction,  end  the  tangent  at  the  midpoint  was 
always  horizontal.  The  application  of  van  Wi jngaardenT s 
theory  gave  only  one  load  deflection  curve  for  non- 
dimensional  plotting,  whereas  Kennedy  obtained  several 
different,  although  similar,  curves  experimentally  by 
employing  varying  arch  geometries.  He  suggested  that 
this  discrepancy  was  due  to  poor  arch  clamping.  The 
experiments  were  conducted  on  arches  fashioned  from 
banding  steel,  from  0.012  inches  to  0.025  inches  in  depth, 
with  a  radius  variation  of  12  inches  to  30  inches,  and 
upon  arches  of  circular  section,  made  from  piano  wire. 

In  a  recent  text  on  flexible  bars,  R.  Frisch-Fay^ 
summarized  the  van  Wijngaarden  treatment,  together  with 
similar  theoretical  work  on  the  buckling  of  a  semi¬ 
circular  ring  by  Biezono  and  Koch,  completed  at  about  the 
same  time  in  1946. 

The  only  other  approach  to  this  problem  was 
that  employed  by  Gjelsvik  and  Bodner^,  who  presented  an 
energy  methods  solution  as  late  as  1962,  which  is  discussed 
at  length  in  a  later  section  of  this  thesis.  Their  theory 
covered  asymmetrical  as  well  as  symmetrical  buckling  modes 
and  related  the  relative  occurrence  of  both  to  a 
dimensionless  parameter  of  the  arch. 
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In  the  same  paper  experimental  evidence  was  also 
presented  which  compared  favourably  with  the  theoretical 
solution.  Load  deflection  curves  were  included  for  some  of 
the  specimens  they  used,  but  the  results  from  the 
experimental  work  do  not  seem  to  have  been  adequately 
evaluated.  The  arch  specimens  were  manufactured  from 
aluminum  alloy  one  inch  wide  and  inch  thick,  and  the 
desired  circular  profile  was  achieved  by  rolling. 

With  no  confirmatory  experimental  consideration 
on  Gjelsvik  and  Bodner’s  work  to  date,  the  author  at  the 
instigation  of  Dr.  J.  S.  Kennedy,  felt  it  worthwhile  to 
examine  thoroughly  and  re-appraise  this  subject.  Gjelsvik 
and  Bodner’s  approach  was  adopted  for  the  theoretical 
analysis,  as  it  had  the  advantage  of  encompassing 
transitional  mode  buckling. 

The  van  Wijngaarden  method,  evolved  before 
digital  computers  were  in  common  use,  has  its  solution 
in  elliptic  integral  form.  This  is  most  difficult  to 
incorporate  in  computer  programmes,  because  of  the 
necessity  to  store  in  the  machine’s  memory  sections  many 
tabulated  values  from  integral  tables. 

The  Mathematics  Department,  at  the  University 
of  Alberta,  advised  that  a  direct  trial  and  error  method 
of  substituting  numerical  values  should  be  employed  in 
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the  attempt  to  satisfy  van  Wi jngaarden1 s  equations. 

Such  a  procedure  had  been  adopted  previously  by  Kennedy, 
and  the  curve  he  obtained  is  reproduced  in  Chapter  4 
for  comparison  with  Gjelsvik  and  Bodnerfs  theory. 

Further  modification  of  this  method,  to  include  the 
possibility  of  an  asymmetrical  buckling  mode  resulted 
in  non-linear  equations,  the  solution  of  which,  thus  far, 
it  has  been  found  impossible  to  determine. 
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CHAPT3R  2 


THEORY 

Since  some  of  the  conclusions  Gjelsvik  and 
Bodner  presented  in  their  treatise  on  snap  buckling  are 
questioned  by  this  author,  in  the  light  of  their  own  theory 
their  theoretical  considerations  are  given  in  this  chapter. 
(A)  General  Formulation 

Consider  a  circular  arch  with  built-in  ends  as 
shown  in  Figure  1.  For  a  thin  curved  beam  the  axial 
strain  is  given  by 

e  =  -  w)  +  _lfdw\2  (1) 

R  doc  2R\doc/ 


where  v  and  w  are  the  tangential  and  radial 
displacements  respectively. 


k,  the  change  in  curvature,  is  given  by 


The  strain  energy  due  to  axial  deformation 


(2) 


is 


U«.  =  b 


-J3 


e  doc 


(3) 


where  Um  has  been  put  in  non-dimensional  form  by 
dividing  by  the  factor  EtbR. 

In  a  similar  non-dimensional  form  the 


strain  energy  due  to  bending  is 
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FIG,  5 


CLAMPED  CIRCULAR  ARCH  UNDER  CONCENTRATED  LOAD 


The 


24  J-£ 

change  in  potential  energy  of  the 
load  which  has  undergone  a  vertical  displacement 
is 

Up  =  -  — L  (w) 
v  EtbR  0 

Hence  the  total  energy  of  the  system 
may  be  expressed  as 


H 


U  +  U  + 
m  p  b 


Now  for  equilibrium 


dH  =  0 
dv 

Since  v  appears  only  in 
strain  energy,  xve  have 


Um  * 

re  9 

h  e  doc 

dH 

=  Mm 

=  Mm  M  = 

Mm 

dv 

dv 

de  d  v 

de 

Um,  the  membrane 


de  doc 
dcx  dv 


=  Mm  M  1 

de  r<oc  (dv/doc) 


— 

rjs  ,  1 

r  e  occ 

de 

h  J 

_doc_ 

LeR  J 

Since  e  will  be  a  symmetrical  function,  the  first 
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bracket  is  not  zero.  The  strain  cannot  be  infinite, 

so  that  the  third  bracket  cannot  be  zero.  Thus  for 

the  equilibrium  condition  to  be  fulfilled  we  must  have 

de  =  0 
doc 

so  that  e  =  C 

where  C  is  a  constant  term. 

Thus  re-arranging  equation  (1) 

2 

dv  =  w  -  1  / dw\  -+  CR 

doc  2R  \doc  j 

Integrating  between  the  limits  <x  = ,  and 
employing  the  boundary  conditions  (v  =  0  at  the  ends 
of  the  arch)  yields 


e  =  C  s  -  1 


2R£ 

Equation  (6)  becomes 


-fVw  -  1  /dw ^  \dpc 

-jo,  \  2R\d<?c  i  / 


(7) 


(a) 


h  =  J 


2J  2 
e  doc  +  t 


2 

k  doc  -  _R 


f 


(w) 

EtbR  oc  =  0 


so  that  substituting  for  e  from  {8)  and  for  k 
from  (2),  H  can  now  be  expressed  in  the  form 


(9) 


H  = 


P 

■fi 


w  -  1 


dw\ 
2R  ^dcc/ 


2i 


d«l2 


/ 


.2 


4- 


24FA 


'  B  /  2  N  2 

*  /dw  \  doc  -  P  (w) 
-JB  Vdoc  j  EtbR 


oc  =  0 


(10) 
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(B)  Determination  of  Buckling  Loads  by  the  Classical 


Eigenvalue  Theory 

The  buckling  load  is  taken  as  the  load  for 
which  geometrically  adjacent  equilibrium  states  can 
exist  in  addition  to  the  original  state  at  the  same  load. 


Consider  small  displacements  about  an  equili 


brium  position  at  the  centre  of  the  arch,  such  that  with 
the  application  of  a  further  small  deflection,  the  arch 
will  change  from  the  unbuckled  to  the  buckled  state. 


Denoting  the  initial  state  prior  to  buckling 


of  the  arch  bv  a  single  superscript,  and  pursuing  a 
similar  line  of  thought  to  that  outlined  in  the  previous 
section,  we  have 


(11) 


Similar  expressions  to  those  obtained  in 


the  previous  section  can  also  be  obtained  for  Hf ,  e ' 
and  k 1 . 

Now  applying  the  double  superscript  to 
denote  the  change  in  these  quantities  during  the 
transition  from  the  unbuckled  to  the  buckled  state, 


o  '  '•  '  ;  r  ■■  '• 

<  :  ~  t,i  "• 

.  n 

'r  f<  t  ; 
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then , 


so  that  the  change  in  total  energy  is 

H"  =  U”  +  U”  +  U" 
m  p  b 

For  equilibrium  in  the  buckled  state 

de**  =  o 
doc 

so  that  upon  integration 

e"  =  C" 

Further,  adapting  equation  (1),  the  change  in  axial 
strain  during  buckling  is 


(12) 


(13) 


e"  =  1 

dv”  - 

w” 

+  -2 

/dw” 

\2  +  2  dw'dw” 

R 

doc 

- 

2R2 

.  vdoc  j 

f  doc  doc  _ 

(14) 


Hence  solving  for  dv” 

doc 


dvn  =  wn  - 

doc  2R 


dw”\2_^_  2  dwTdw” 

da  /  doc  dec 


(15) 


Integration  of  the  above  expression,  and 
a  consideration  of  the  boundary  conditions  yields 


'  5 


- 
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' 
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c  ■ 


li 


C"  _  _ 


2Rfi 


w 


it  _ 


2R 


/dw"\  2  dwf  dw" 

IdoC  /  doc  doc 


►doc 


The  difference  in  the  curvature  change  is 
k" 


1  dw" 

R2^ 


Thus  substituting  (l6)and(17)  into  equation  (12) 
the  total  energy  change  is 


H"  =  1 


4R2£ 


Ww"  -_1 

'A  2R 


/dw"\2  |  2dwT  dw" 

Vdoc  /  doc  doc  . 


doc 


1 2 


(16) 


(17) 


(18) 


'PU 

-J3l  2R fi  J 


'?  ( w"  -JL 
-p  \  2R 


/dw"\  +  2dw*  dw" 
\doc  )  doc  doc 


|docjdoc 


+ 


2 


24R4 


/dw"  \  *  doc  ^  t 

fW2) 


12R2J 


P 

f 


kT  dw"  doc-  P 
2 


(w") 


cc  = 


s  0 


doc' 


EtbR 


Neglecting  terms  in  w  which  are  of  greater 
order  than  second  and  noting  since  the  initial  state  is 
an  equilibrium  state  all  linear  terms  in  w  must  cancel, 
equation (1$)  reduces  to 


H"  - 


4R2j3 


£/w"  -  1  dwT  dw"\doc 
-J3  \  R  da  doc  / 


4- 


2R 


1  (~j3  /dw"\  da  t 
-j3  \doc  / 


2  rn  .  2  v  2 


P  /dw"\  aoc 


(19) 


24R4J-#  Vdoc2/ 


A  general  deformed  shape  may  be  composed  of 
a  symmetrical  function  of cc ,  w^ t  and  an  asymmetrical 
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function  ofoc,  W2.  So  we  may  write 


wn  =  A”  vrj  +  B”  w2 

The  functions  are  thus  orthogonal  in  the  range 


-  J3  ^  oq  <  p 


Hence  (19)  becomes 


H"  =  1 

TsTp 


P 
f  L 


Anwi  4-  BnW2 


1  dw’  f  AT,dwp 

R  doc  \  doc 


+ 


(20) 


Because  of  the  orthogonality  properties  Hn 


can  again  be  reduced  to 


H"  ,  At,2f  i 

fifal  ~ 

ldwT 

dwp^  doc 

„  V. 

"A 

Rdoc 

dec  / 

-i2 


+ 

+ 


e  T 

2R2 

B" 


t2  \»(M 
2*1*  J  fi  [doc2) 


For  eouilibrium  in  the  buckled  state 

3H"  =  0  3H"  =  0 

d  A"  5Bff 


(21) 


Hence  each  of  the  two  expressions  in  parenthesis  in 
(21)  must  be  zero. 
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From  the  second  equilibrium  expression 


rp 

,  2 
rdw2 

i2 

doc 

e '  =  -  X?- 
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or,  if  written  in  terms  of  x  =  oc 
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(22) 
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(23) 


Examining  the  form  of  (23)  above,  it  is  seen 
that  the  force  is  equal  to  the  second  mode  Euler  load 
of  a  clamped  column  of  the  same  length  as  the  arch. 5 

Substituting  from  (3)  for  e’,  the  deflection 
at  buckling,  which  will  be  when  the  transitional  or 
asymmetrical  buckling  load  is  reached,  is  obtained  from 
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Introducing  non-dimensional  amplitude 


an(^  the  parameter  X  = 


the  above  expres¬ 


sion  becomes 
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-I  \dx  J 
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-I  3X2- 
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m2)  dx 


dw2\  2 
dx  J 
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(25) 


dx 


from  which  a  can  be  evaluated  provided  the  displacements 
W1  and  W2  can  be  determined.  This  value  of  a  is  actually 
aT}  the  displacement  at  which  transitional  buckling 


occurs* 


The  first  equilibrium  equation,  obtained  from 


(20)  by  setting 
gives  the  following  relationship 
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(26) 


which  gives  an  expression  eventually  for  the  deflection 
of  the  arch  ay  at  the  upper,  or  symmetrical,  buckling  load. 

This  line  of  thought  was  elaborated  further 
by  Gpelsvik  and  Bodner  and  the  theoretical  results  they 
obtained  were  compared  with  the  results  of  another  approxi¬ 
mate  method  (which  is  contained  in  the  same  paper). 

The  algebraic  expression  for  the  various  loads 
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and  deflections  agreed  exactly,  Using  these  results  they 
derived  a  series  of  relationships  which  for  varying  X 
would  enable  one  to  predict  whether  transitional  buckling 
would  occur  at  all,  and  if  it  did  whether  the  transitional 
mode  would  occur  before  or  after  the  occurrence  of 
symmetrical  buckling.  The  various  possibilities  for  arch 
buckling  are  described  diagrammatically  in  Figure  2. 
Gjelsvik  and  Bodnerrs  contention  was  that  for  sufficiently 
low  values  of  X,  no  buckling  as  such  xrould.  occur  at  all 
for  increasing  central  deflection,  and  the  resulting  load 
deflection  curve  would  have  no  unstable  portion.  As  the 
value  of  X  was  increased,  symmetrical  buckling  only  would 
occur.  Further  increases  in  X  would  initially  cause  a 
transitional  mode  of  asymmetry  to  occur  after  symmetrical 
buckling  as  such  had  occurred,  and  eventually  conditions 
would  be  attained  whereby  the  transitional  buckling  mode 
would  occur  first.  In  this  case  of  course  Py  would  never 
be  reached  for  the  arch  would  have  buckled,  and  the  subse¬ 
quent  values  of  P  would  decrease. 

By  the  nature  of  the  equation  for  the  transi¬ 
tional  mode,  it  was  shown  by  Gjelsvik  and  Bodner  that  a 
straight  line  would  always  indicate  such  a  mode.  For 
exceptionally  large  values  of  X  transitional  mode  buckling 
could  happen  shortly  after  the  load  had  risen  to  the 
intermediate  buckling  load  value;  the  intermediate 
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SYMMETRICAL  3UCKLING 
WITH  TRANSITIONAL  MODE 
PRESENT  IN  UNSTABLE  REGION 


TRANSITIONAL  ASYMMETRICAL 
BUCKLING 


FIG.  2 


VARIATION  OF  BUCKLING  MODES 
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buckling  load  was  shown  to  be  a  lower  bound  for  snap 
buckling,  which  is  buckling  involving  a  large  deforms 
tion  before  another  equilibrium  position  is  attained. 


Unfortunately  these  methods  and  the  result¬ 


ing  criteria  for  the  form  of  buckling  are  difficult  to 
check  numerically.  The  author  could  not  find  an 
asymmetrical  function,  W£,  which  would  satisfy  the 
boundary  conditions  required  for  the  built-in  ends,  and 
no  such  function  was  suggested  by  Gjelsvik  and  Bodner. 


Therefore  only  sufficient  of  the  above  theory 


has  been  developed  that  is  necessary  for  the  following 
series  approach. 

(C)  Series  Solution  of  the  Non-Linear  Deformations  of 
the  Clamped  Arch. 

Consider  initially  a  symmetrically  deflected 
shape,  and  let  this  shape  be  described  by  the  function 


r 


(27) 


where  wn  is  a  symmetrical  function,  and  An  is  the 
amplitude  of  this  function. 

Substituting  the  deflected  shape  into 
eauation  (10)  the  total  energy  may  be  expressed  as 
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(23) 


where  an  s  An 

and  H  *  =  HX 

W 

F 

wn  should  be  selected 

so  that  dwn  and 

,2 

dwn  form 

dx 

dx2 

(29) 


orthogonal  sets  such  that 
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-I  dx  dx 
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Then  H  can  be  reduced  to 
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The  equilibrium  condition  is 


dH  s  0 
dan 


(31) 


which  then  renders  r  equations  of  the  form 
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After  Gjelsvik  and  Bodner  the  following 
notation  is  introduced 


'In  = 


wndx 


'2n  = 


'3n  = 


2 

*  dx 

-iVdx2; 


J-|\dx 
C  =  (wn) 


^n\2dx 


An  '  n '  x  =  0 
Thus  modifying  eouation  (32) 


X 
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r  r  r 

^~ancln  “  ?  ^anc2n 


Gln  -  anC2n 
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+  ib  an°3n " F  °4n  =  0 


The  set  of  equations  (34)  are  the 
equilibrium  equations  for  symmetrical  deflections 
of  the  arch.  These  equations  are  later  transformed 
in  order  that  P  can  be  plotted  against  a. 


(32) 


(33) 


(34) 


' 


Now  applying  the  condition  obtained  by  the 
classical  eigenvalue  considerations  for  asymmetri cal 


buckling  to  occur,  that  is,  expression 
keeping 


w  - 


V'n<x> 


I 


{21+)  ,  but 


then 


where 


K  = 


% 

dx 

3T 

dx 


(35) 


(36) 


The  constant  K  is  the  coefficient  in  the 
second  mode  buckling  equation  for  a  clamped  column^ 

so  that 

K  =  2 . 0l+26rC 

Here  wfi  is  the  asymmetrical  buckled  shape. 

Now  substituting  the  conditions  defined  in 

(33)  the  criterion  for  asymmetrical  buckling  becomes 

r  r 

y"anCln  "  i  Y  anc2n  =  _ K 

.  .  6F 

Hence  to  obtain  the  asymmetrical  buckling 
load  Pip  the  right  hand  side  of  equation  (37)  is 


(37) 


•  ■  • 

♦ 

- 

inserted  into  (34)  yielding 


Zi 


Yfx[Cln  ‘  SnC2n  +  jfj  a"C3n  -  PTCAn  =  0 


(38) 


And  so  anf(  the  non-dimensional  deflection 


at  P-pjis 


anT  =  12  I  PT  C4n  -  KCln 


(39) 


c3n  -  KC2n 

_  * 

If  P-p  is  known,  then  this  expression  may 
be  applied  directly  to  evaluate  an«p. 

Reworking  (37),  by  replacing  an  with 
expression  (39)  above,  the  criterion  for  asymmetrical 
buckling  becomes 
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Introducing  the  following  notation 


■ 


22 


y  rcln  c3n  ~ 

4  LC3n  -  KC2n)2- 


D3  *  f  fC?n  (2C^n  -  KC2n)] 

XL  (C3n  -  KC2n)2  J 

the  asymmetrical  buckling  load  is  given  in  the 
following  equation 


(12  X  P*)  2Pp  -  24  X  P^D2  +  JC  -f-  KDo  =  0 


3X 


2 


so  that 


(41) 


(42) 


Therefore  if  a  satisfactory  function  wn 
% 

can  be  found  P^  may  be  calculated  after  initially 

determining  the  various  G  and  D  values,  and  its 

corresponding  an<p  value  may  be  evaluated  from 

(39).  To  determine  P^  it  is  necessary  to  apply  equations 

(3X),  and  by  plotting  P  against  a,  Py  will  be  given  at 

the  maximum  value,  which  terminates  the  first  stable 

section. 

These  steps  will  be  elaborated  on  in  the 


next  section. 
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( D )  Application  of  a  General  Displacement  Function  to 

the  Series  Solution. 

Consider  the  following  function 
wn  -  i  (1+cos  nrrx)  ,  n  odd 
wn  =  h  (1  -  cos  rmx)  ,  n  even 


F  I  G0  3  SKETCH  OF  FUNCTION  [43]  FOR 
INITIAL  VALUES 


This  function  satisfies  the  boundary  conditions 
imposed  by  clamping  the  arch  at  its  ends  and  also  it  meets 
the  conditions  of  orthogonality  required  by  equation  (30). 
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-  all  n, 


-odd  n, 

-  even  n, 


Determining  now  the  various  C  and  D  values, 
cln  =  1 
c2n  =  n2TT2/4 
c3n  -  ni*TtV4 
c4n  «  1 


'4  n 
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The  ensuing  D  terms  are  thus 
r  r 

D1  =  I  nVlfr  -  2.0426 ) 2 


,  odd  n 


(44) 
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-  ,  all  n 


Hence  P*j>  may  be  calculated  for  solutions 


containing  a  varying  number  of  terms.  In  their  paper 
G.jelsvik  and  Bodner  stated  that  although  a  substantial 
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difference  is  obtained  between  the  one  and  the  two  term 
solution,  very  little  deviation  is  recorded  between  the 
latter  and  a  six  term  solution. 

However  the  application  of  a  series  with  more 
than  two  terms  to  evaluate  the  upper  buckling  load  is 
an  extremely  arduous  undertaking  and  subseouently  the 
two  term  solution  will  be  derived.  This  is  later  shown 
to  be  of  sufficient  adequacy. 

Taking  r  =  2,  from  equation  (42) 


P 


T  = 


TT 

12X 


1  +-  (0.2245  -  0.1&47  ttM 

X2 


(46) 


To  obtain  the  load  deflection  curve  and  the 
upper  buckling  load,  equations  (34)  are  now  developed. 
For  any  number  of  terms,  the  expression  below  follows 
immediately 
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For  r  a  2,  equation  (47)  reduces  to 


a2  =  1 


-  (iriu) 


j  m  * p 


(48) 


2  Y 

where  d  =  tt  ai  and  u  = 

TT^ 

Thus  in  the  first  of  the  equilibrium 
eauations  (34),  a2  can  ^e  replaced  by  (4$),  which 
results  in  an  expression  relating  F  and  a.  In  a 
two  term  series  solution,  it  should  be  noted  that  a, 
the  non-dimensional  central  deflection  of  the  arch, 
is  equal  to  a^.  This  is  adequately  clarified  by 
Figure  3. 

Developing  the  first  equilibrium  equation 


2  2  2  4  o  ,  -| 
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(49) 


Substituting  for  a2  from  (4$) 

(d  -  d2)  (h2  t  P*)2  -  (hi  -  P*)  (h2  +  P*)  -  |  (hx  -  F^)2 
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(50) 
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The  load  deflection  relation  is  thus 
determined,  although  in  a  particularly  awkward 
fashion,  by  (50)  above.  For  simplification  consider 
the  new  variables 


a  =  d  -  K 

Rewriting  (50),  we  have 


p  =  4-SuP  -  4 


Ko  +  Knp  +  K?p^  -  4p3  =  0 


whe  re 


-3 


K0  =  a 


Ki  = 


36  +  3u2  (9a  -  116) 


-12  +  3u2  (6^  -  120) 


* 


(51) 


(52) 


(53) 


K2  =  a  “20  -+-  3u^  (a^  -  20) 

With  the  aid  of  (52)  values  of  p,  and  thus 
P^from  (51),  may  be  obtained  for  varying  d  values  for 
any  given  X. 

G.jelsvik  and  Bodner  did  not  state  how  they 
obtained  successive  values  from  (52),  which  is  a  lengthy 
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e 


process  by  manual  calculation. 

A  computer  programme  was  devised  (Appendix  B) 

* 

which  rendered  the  successive  values  of  P  against  ap 
for  a  wide  range  of  X  values.  The  X  values  selected 
included  those  in  the  range  which  Gjelsvik  and  Bodner 
examined  experimentally  and  theoretically,  through  to 
quite  large  values. 

The  load  deflection  curves  from  the  results 
furnished  by  the  University  of  Alberta  IBM  1620  computer 
are  presented  in  Chapter  A,  where  they  may  be  simul¬ 
taneously  compared  with  the  experimental  plots. 

Also  depicted  in  the  same  section  is  the 
theoretical  curve  of  load  against  deflection  by  applying 
van  Wi(jngaardenT s  theory,  which  was  evaluated  by 
J.  S.  Kennedy  at  an  earlier  date. 
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CHAPTER  3 


APPARATUS  AND  EXPERIMENTAL  PROCEDURE 
(A)  Aparatus 

The  experimental  equipment  used  is  shown  in 

Figure  4. 

The  arch  bed  was  a  length  of  65S-T6  6nx2^" 

aluminum  channel,  weighing  4.07  pounds  per  foot;  the 
top  face  and  the  tips  of  the  channel  side  walls  were 
milled  in  the  department  workshop  to  ensure  that  these 
surfaces  would  be  smooth  and  parallel  to  each  other. 

The  top  face  of  the  channel  was  then  machined 
as  in  Figure  5.  The  centre  slot  would  thus  enable 
sufficiently  large  deformations  to  be  applied  to  the 
arch,  should  it  prove  necessary  during  a  test  for  the 
arch  centre  to  be  deflected  to  a  very  low  position. 

The  parallel  rows  of  holes  were  drilled  so  that  the  arch 
clamp  base  would  slot  in  via  guide  pegs,  and  would  thus 
be  aligned  correctly  with  no  further  effort.  The  holes 
tapped  at  regular  intervals  along  the  channel  centre¬ 
line  were  so  placed  to  mate  with  the  fixing  bolt  of 
each  clamp  base. 

The  arch  clamp  surfaces  were  milled  so  that 
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FIG.  4  GENERAL  VIEW  OF  APPARATUS 
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PLAN  VIEW  OF  ARCH  BED 


when  the  clamps  were  set  in  position  on  the  arch  bed, 

the  angle  of  inclination  of  the  arch  built-in  end 

0 

would  be  at  30  to  the  horizontal.  Details  of  the  arch 
clamp  structure  are  depicted  in  Figure  6. 

The  arch  bed  was  fixed  to  the  scale  pan  of  a 
Troemner  field  scale,  which  in  turn  had  been  securely 
located  on  a  firm  timber  bench.  The  scale  pan  had  been 
skimmed  on  a  lathe  to  ensure  that  the  arch  bed  would 
rest  squarely  upon  it.  Careful  adjustment  of  the  bed 
scale  system  now  followed  and  it  was  checked  with  a 
sensitive  spirit  level  to  ensure  that  the  upper  surface 
of  the  channel  lay  in  the  horizontal  plane.  The  eventual 
degree  of  horizontal  register  compared  favourably  with  the 
Department  milling  machine  bed. 

Subsequent  arch  deflections  were  procured  by 

a  hook  crau^e,  which  was  suooorted  by  a  tubular  steel  frame. 

The  vertical  alignment  of  the  gauge  measure  was  set  with  the 
aid  of  a  field  transit.  Two  interchangeable  loading  heads 
were  constructed  to  fit  squarely  on  the  bottom  of  the  gauge 
measure;  a  loading  head  which  would  permit  any  of  the  ensuing 
test  arches  to  deform  only  in  a  symmetrical  mode,  as  in  Figure 
7  and  3,  and  an  asymmetrical  loading  head  which  would  not 
restrict  the  manner  of  the  deformations.  Figures  9  to  11 
illustrate  the  latter. 
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FIG.  6  CLOSE-UP  OF  ARCH  CLAMP 
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FIG.  7  SYMMETRIC  HEAD,  NO  LOAD 
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FIG.  8  SYMMETRIC  HEAD,  ARCH  BUCKLED  SYMMETRICALLY 
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FIG.  9  ARCH  PRIOR  TO  BUCKLING  RETAINING  SYMMETRICAL  MODE 


FIG.  10  ARCH  IN  UNSTABLE  REGION,  EXHIBITING  PRONOUNCED  ASYMMETRY 
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ARCH  DEFLECTED  TO  SECOND  STABLE  REGION,  SHOWING  RETURN  TO  SYMMETRY 


Care  was  taken  to  place  the  scale-arch  bed 
system  and  the  gauge  support  frame  relative  to  each 
other,  so  that  final  adjustment  would  leave  the 
loading  head  accurately  lined  up  with  the  arch  bed 
centre-lines. 

Prior  to  its  installation,  the  Troemner 
field  scale  had  been  investigated  for  accuracy,  using 
weights  checked  on  a  Mettler  comparater  scale.  The 
scale  cumulative  error  was  found  to  be  a  maximum  of 
±  5  grams, 

(B)  Test  Specimens 

All  the  arch  test-pieces  were  cut  from 
banding  steel  kindly  supplied  by  Acme  Steel,  The 
circumferential  arch  lengths  for  the  varying  clamp 
settings  were  scribed  by  a  milling  machine,  in  each 
case  leaving  sufficient  overlap  to  be  accomodated 
comfortably  in  the  clamps.  In  addition  centre-lines 
were  scribed  on  each  piece  to  correspond  to  similar 
marks  on  the  clamps  and  thereby  assisted  in  alignment. 
The  painted  black  finish  of  the  banding  steel  lent 
itself  readily  to  this  process.  Each  length  was 

diligently  inspected  for  imperfections  and  kinks. 

Seven  arches  were  considered  experimentally; 


■  ••  o.':-  rjf. 

*  -!• 

'  • r  • , 

■ 

.  ■ 

...  fr. 

. 

r  —  r-  't  '■n  ■  < 

.  (  ' 

«*  r  •  r. 

■  -  ps  e 

f,  • . ■  $'  r  riorA  •  t  j.r.i 

I  :  18 

' 

. 

u  u'*"  rv-J  -  V  >  3twt0  •  t*l 

'  :•  ■  tSlbba.  .  f  I  ■ 

’  ,f'*  p  r>  n  .-n-  7  -R  b,  d  ;-v>c  :  t«w 

j  .t;  »n  or:  f 

.  ' :  : '  0  1  i  ‘  ‘ 

.  ■  **  f-,  • 

.5rr  mlT'3  !X<  bei:’hiwo  n#w  a  >v«?c 


dimensions  are  given  in  Table  1  and ,  for  the  sake  of 
convenience,  each  specimen  has  been  given  an  identity- 
number  . 


TABLE  1.  TEST  SPECIMENS 


IDENTITY 

NO. 

SPAN 

DEPTH 

WIDTH 

X  =  £2R/t 

SI 

15” 

0.0152” 

0.75" 

270.5 

S2 

20” 

0.0149” 

0.75” 

369.0 

S3 

25” 

0.014L” 

0.75" 

475.9 

S4 

30” 

0.0155” 

0.75” 

530.6 

S5 

20” 

0.0197” 

0.625” 

273.3 

S6 

20" 

0.0210” 

0.5  " 

261.1 

S7 

30” 

0.0127" 

0.5  " 

647.6 

The 

X  values  of  these 

arches  were 

all  quite 

large,  compared  to  the  range  Gjelsvik  and  Bodner 
investigated,  which  varied  from  X  -  3.69  to  X  =  16.25. 
Although  it  was  felt  desirable  to  include  specimens 
which  did  have  smaller  X  values,  the  correspondingly 
higher  buckling  loads  which  would  result  were  outside 
the  scope  of  the  Troemner  scale.  But  had  a  large 
scale  been  employed  initially,  then  a  considerable 
loss  of  accuracy  would  have  been  sustained  in  the  low 
load  readings. 
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(C )  Experimental  Procedure 

With  a  test  piece  firmly  located  by  the  clamps 
the  jockey  weight  of  the  scale  was  adjusted  so  that  the 
scale  reading  was  zero  under  the  action  of  the  undeformed 
arch  and  the  arch  bed.  The  hook  gauge  measure  was  now 

lowered  until  the  loading  head  just  touched  the  arch. 

In  the  case  of  the  symmetrical  head,  the  head  was  securely 
fixed  to  the  arch  centre. 

The  measure  was  then  deflected  incrementally, 
and  for  each  setting  the  balancing  weight  registered  on 
the  scale  arm  was  noted;  hence  with  the  scale  adjusted  so 
that  the  scale  pan  was  again  held  at  the  original  height, 
a  situation  was  obtained  whereby  a  known  central  load, 
induced  the  given  arch  to  deflect  a  measured  distance. 

For  each  specimen  the  asymmetrical  loading  head 
was  set  onto  the  hook  gauge  measure  immediately  after  a 
test  with  the  restrictive  head,  to  facilitate  direct 
comparisons.  The  symmetrical  head  was  replaced  afterwards 
to  check  the  scale  zero  reading  for  consistency.  With  the 
advent  of  any  such  deviation,  the  recently  acquired 
readings  were  rejected. 

Precautions  were  taken  to  load  the  arch  with  the 
free  head  very  gently  in  the  region  of  the  upper  buckling 
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load  as  sudden  deflections  would  possibly  cause  the  arch 
to  buckle  transitionally  under  the  influence  of  the 
sudden  motion,  rather  than  the  static  load.  Whenever 
the  arch  was  deformed  asymmetrically,  the  point  of 
application  of  the  load  was  moved  until  it  was  again  at 
the  mid-length  of  the  arch.  However,  lateral  movement  of 
the  head's  roller  registered  little,  if  any,  difference 
with  the  scale  and  so  the  roller  was  kept  slightly  to 
the  inner  side  of  the  centre  scribe  line.  This  minor 
precaution  prevented  inconvenient  snap-throughs  of  the 
asymmetry  to  the  opposite  side  of  the  arch. 

Results  recorded  in  this  manner  were  repeatable 
to  within  the  scale  error.  Difficulty  was  experienced 
with  the  shortest  arch  length,  where  it  was  found  that 
repeated  discrepancies  were  caused  by  yielding  at  the 
ends.  By  avoiding  any  further  deflection  just  after 
the  lower  buckling  load  was  attained,  this  problem  was 


overcome 
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CHAPTER  4 


RESULTS  AND  DISCUSSION 


( A )  Symmetrical  Buckling 

Dimensional  analysis  of  the  problem  revealed  a 
function  containing  the  following  possible  dimensionless 

PrZ-2 

groups,  —  , ,  where  z  is  any  number,  A,  b,  R  ft  and 

Ebt  R  t  t,r 

Poisson* s  ratio.  By  the  use  of  only  one  material  for  the 

test-pieces,  Poisson* s  ratio  was  kept  at  a  constant  value; 

all  the  experiments  were  conducted  for  a  particular  angle 

subtended  at  the  centre  of  the  circular  arch,  which  fixed  J3. 

This  left  the  relationship 


Neither  van  Wijngaarden  nor  Gjelsvik  and  Bodner  indicated 
that  the  variation  of  the  breadth  to  thickness  ratio  of 
the  arch  would  have  an  effect  on  the  resulting  load  deflec¬ 
tion  curves.  To  examine  any  possible  influence  of  this 
dimensionless  parameter,  the  load  deflection  curves  for 
three  specimens  SI,  S5,  and  S6,  which  had  an  approximately 

common  R  ratio,  have  been  plotted  in  Figure  19.  Van  Wijn- 
t 

gaarden  coordinates  were  chosen  for  later  convenience, 


otherwise  the  effect  could  have  been  examined  as  easily  on 
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a  G jelsvik-Bodner  basis  whereby  X  =  was  held  constant. 
_  t 

From  thL  s  diagram  it  may  be  observed  that  any  variation 

in  the  parameter  b,  over  the  range  involved  in  the  tests, 

t 

seems  to  have  no  predictable  influence  upon  the  relative 
position,  or  shape,  of  the  curve. 

Theoretical  load  deflection  curves  for  the 
four  pieces  SI  -  S 4  of  varying  X  values,  obtained  by 
application  of  the  G jelsvik-Bodner  theory,  have  been  plot¬ 
ted  in  Figure  12.  On  this  non-dimensional  basis  it  may  be 
seen  that  the  upper  and  lower  buckling  loads  occurred  at 
consistent  deflections,  a  values  of  0.155  and  0.65 
respectively.  As  the  value  of  X  was  increased,  the  buck¬ 
ling  loads  for  each  X  decreased;  this  effect  is  further 
illustrated  in  Figure  17. 

The  experimental  results  for  these  X  values  are 
depicted  in  Figures  13  -  16,  plotted  in  G jelsvik-Bodner 
coordinates,  showing  values  obtained  with  the  symmetrical 
and  asymmetrical  heads.  The  latter  are  discussed  in 
Section  (B)  of  this  chapter.  The  theoretical  and  experi¬ 
mental  values  for  a^  and  a^  agreed  quite  closely  for  the 
symmetrical  case.  However,  there  was  a  consistent  dis¬ 
crepancy  evident  whereby  each  experimental  curve  was  lower 
than  the  theoretical  counterpart,  that  is  for  any  given 
deflection  the  experimental  load  value  was  less  than  the 
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theoretical. 

Developement  of  the  van  Wijngaarden  theory- 
yields  only  one  load  deflection  curve,  and  by  transforming 
the  experimental  load  displacement  values  for  specimens 
SI  -  S4  to  van  Wijngaarden  coordinates,  the  curves  given 
in  Figure  IS  were  acquired.  From  this  figure  it  may  be 
seen  that  as  the  radius  to  thickness  ratio  of  the  arch 
was  increased,  the  curves  were  correspondingly  lowered. 

In  each  case  the  experimental  curve  was  lower  than  the 
theoretical  van  Wijngaarden  curve.  Thus  it  appeared  that 
possibly  van  Wijngaarden,  in  his  solution,  unlike  Gjelsvik 
and  Bodner,  did  not  make  allowance  for  the  parameter  R. 

Investigating  further  this  seeming  lack  of  a 
parameter,  the  theoretical  results  obtained  for  the 
specimens  SI  -  S4  were  transferred  to  van  Wijngaarden 
coordinates,  and  a  single  load  deflection  curve  was 
obtained.  This  is  shown  in  Figure  22,  where  it  may  be 
compared  with  the  curve  given  by  van  Wijngaarden* s  theory. 
The  only  deviation  at  all  between  these  two  curves  is  a 
difference  of  2  in  P,  at  the  lower  buckling  load,  and 
otherwise  they  concur  most  favourably. 

Kennedy  calculated  the  following  values: 

Py  =  2&.7  P^  =  10.2  occurring  at 
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L  =  0.170 


=  0.041 
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The  transformed  Gjelsvik-Bodner  values  yielded 
Py  =  2£.73  Pl  "  3.96  occurring  at 

& u  =  0.0417  &L  =  0.180 

Following  on  from  this  comparative  analysis, 
the  discrepancies  between  the  theoretical  and  experimental 
results  were  still  unaccounted  for  satisfactorily.  Since 
the  two  theories  were  in  agreement  it  was  concluded  that 
the  experimental  results  must  be  in  error  and  it  was 
suspected  that  the  weight  of  the  arch,  neglected  in  both 
theories,  was  the  c ause  of  the  consistently  low  results. 
Subsequent  weighing  of  a  length  of  banding  steel, 
lf  x  0.75”  x  0.015”,  lent  support  to  this  line  of  thought, 
as  its  weight  was  17  grams.  On  this  basis  the  arch  weight 
for  pieces  SI  -  S4  have  been  calculated  and  are  given  in 
Table  2. 

TABLE  2.  ARCH  WEIGHT-  BUCKLING  LOAD  RATIOS 


SPECIMEN 

X 

BUCKLING 

LOAD  (gms. ) 

ARCH 

LENGTH 

WEIGHT 
(gms. ) 

WEIGHT/Py 

SI 

270.5 

354 

15.703” 

22.2 

.063 

S2 

363.0 

17# 

20.944" 

29.6 

.366 

S3 

475.9 

96 

26.130” 

37.1 

vO 

•CO 

. 

S4 

530.6 

73 

31.416" 

44.5 

.570 

The  weight  to  buckling  load  ratio  increased  considerably  as 
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the  X  value  increased,  which  would  explain  why  the  dis¬ 
crepancy  between  theoretical  and  experimental  results 
became  greater  for  the  larger  X  values. 

From  Figure  16,  it  appeared  that  for  test 
specimens  with  X  greater  than  550,  a  portion  of  the  curve 
would  be  below  the  deflection  axis,  that  is,  the  lower 
buckling  load  would  be  negative.  Test  piece  S7  was  there 
fore  made  up,  and  it  was  discovered  that  this  specimen 
would  indeed  remain  freely  in  the  buckled  state  as 
surmised;  Figure  2L  confirms  this.  This  contradicts 
Gjelsvik  and  Bodnerfs  general  statement  "The  lower 
buckling  load  of  the  clamped  arch  is  greater  than  zero 
and  therefore  the  arch  cannot  remain  in  the  buckled  state 
when  the  load  is  removed."  The  actual  curve  obtained  is 
presented  in  Figure  25,  although  it  should  be  emphasized 
that  it  is  of  dubious  accuracy,  due  to  a  combination  of 
the  poor  sensitivity  of  the  Troemner  scale  and  the  small 
loads  being  measured  (0  to  30  grams). 

Nonetheless,  the  characteristic  difference 
resulting  from  the  use  of  the  two  different  loading  heads 
is  still  preserved.  Such  traits  are  reviewed  in  the  next 
section. 

Gjelsvik  and  Bodner  declared  that  for  a  very 
low  X  ratio,  no  buckling  at  all  would  be  experienced  by 
the  arch,  or,  that  at  no  time  would  there  be  a  region  of 
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instability.  Accordingly  P]j  and  were  calculated  for  a 
wide  range  of  X  values  and  are  stated  in  van  Wijngaarden 
non-dimensional  terms  in  Table  3. 

TABLE  3.  BUCKLING  LOADS  BY  GJSLSVIK-BOBNSR  THEORY 


TRANSFORMED  TO  VAN  WIJNGAARDEN  COORDINATES 


X 

p* 

PU 

pu 

p! 

PL 

3.6° 

0.231* 

19.6 

.2140 

18.1 

5.00 

.1943 

22.2 

.1366 

15.9 

8.01 

.  1402 

25.9 

.06516 

12.0 

11.62 

.1024 

27.3 

.03916 

10.4 

14.50 

.08357 

27.8 

.02Q89 

9.93 

18.00 

.06808 

28.1 

.02330 

9.61 

25.00 

.04954 

28.4 

.01626 

9.32 

30.00 

. 04144 

28.5 

.01339 

9.21 

35.00 

.03560 

28.6 

.01140 

9.15 

40.00 

.03120 

28.60 

.009924 

9.10 

45.00 

.02776 

28.63 

.008993 

9.07 

50.00 

.02500 

28.65 

.007896 

9.05 

100.00 

.01253 

28.71 

.003919 

8.98 

150.00 

.008359 

28.93 

.002609 

8.99 

270.50 

.004635 

28.74 

.001445 

8.96 

320.00 

.003918 

28.73 

.001221 

8.96 

368.00 

.003407 

28.73 

.001062 

8.96 

420.00 

.002985 

28.73 

.000930 

8.97 

475.90 

.002634 

28.73 

.000821 

8.95 

530.60 

.002363 

28.74 

.000736 

8.95 

P  =  12X  P 

where  J3  =  0. 

5236 

J3 

Inspection  of 

this  table, 

and  of  Figure 

23, 

indicates 

that  one  theore 

tical  curve 

can  indeed  be 

fore- 

cast  for 

values  of  X  greater  than  100 

However,  for  lower 

‘  ic  b 


. 


■ 


- 


— r 

<* 

'  1  1  „ 

* 

r 

‘  " 

*» 

, 

.  -  r 

*r 

-» 

.  :• 

' 

T» 

A 

r* 

# 

. 

]•  „ 

1. 

. 

'  :* 

. 

» 

* 

„ 

f  r  r; 

• 

« 

• 

,  r  ,  ■  ,  •  o 

'  1 

m 

* 

•  V'  ij 

. 

1  , 

" . 

. 

. 

*1 

<  . 

• 

• 

. 

• 

. 

. 

* 

. 

. 

•  r 

■* 

« 

. 

C . 

»> 

o;  .  1 

1  ^ 

. 

* 

• 

» 

.  ri 

* 

• 

o 

* 

« 

. 

* 

.  ■ 

. 

. 

. 

• 

• 

. 

. 

t 


t 


« 


values  of  X,  particularly  values  less  than  20,  this  is  not 
so  and  van  Wi jngaarden ’ s  theory  is  not  adequate.  Figure  23 
infers  that  buckling  is  first  encountered  at  an  X  ratio  of 
approximately  3*5,  whence  with  subsequent  increases  in  this 
parameter,  Py  and  diverge  rapidly  to  approach  asympto¬ 
tically  the  values  23.73  and  3.96  respectively. 

At  this  stage  it  should  be  pointed  out  that,  in 

their  paper,  Gjelsvik  and  Bodner  stated  emphatically  that 

around  X  =  10,  the  asymmetrical  mode  became  dominant,  and 

*- 

they  alluded  that  Py  values  became  very  large  after  this 
criterion.  The  author  does  not  believe  this,  and  this 
point  is  discussed  in  the  ensuing  section. 

(B)  Asymmetrical  Buckling 

The  experimental  transitional  modes  recorded  are 
presented  in  the  diagrams  which  contain  the  experimental 
symmetrical  plots,  discussed  previously;  namely  Figures 
13  -  16 ,  20,  21  and  25.  It  was  found  that  employing  the 
non-restrictive  loading  head,  the  arches  initially  displayed 
symmetrical  deformations  until  a  deflection  ay  had  been 
applied  to  the  arch  centre.  At  this  stage  the  load  was 
slightly  less  than  the  upper  buckling  load.  Snap  buckling 
into  the  transitional  mode  occurred  at  about  this  deflection 
and  the  arch  profile  assumed  a  noticeable  degree  of  asymmetry. 
Further  vertical  deflections  gave  a  load  deflection  relation 
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which  in  each  case,  fitted  approximately  to  a  straight  line 
when  plotted.  This  line  passed  through  a  point  near  to  the 
intersection  of  the  intermediate  buckling  load  line  and  the 
symmetrical  load  deflection  curve,  as  predicted  by  theor¬ 
etical  considerations.  For  high  loads  the  points  actually 
lay  on  a  curve  which  did  not  merge  into  a  definite  straight 
line  until  shortly  before  the  aforementioned  intersection. 
Even  so,  the  experimental  occurrence  of  secondary  buckling 
is  clearly  noticeable  on  the  graphs.  With  the  exception  of 
the  S7  plot,  which  was  not  constructed  for  the  negative 
lower  buckling  load  region,  each  asymmetry  plot  rejoined 
the  symmetrical  curve  slightly  before  the  minimum  part  of 
the  curve;  by  this  time,  in  actual  practice,  the  arch  had 
regained  its  symmetrical  shape,  and  further  deflections 
yielded  both  the  arch  shape  and  load  deflection  values 
which  were  the  same  as  those  obtained  with  the  restrictive 
loading  head. 

The  degree  of  variation  experienced  at  the  upper 
load  portion  of  the  curves  fell  within  experimental  error 
allowances.  But  each  arch  tested  had  consistently  slightly 
lower  load  readings  in  this  region,  under  the  action  of  the 
asymmetrical  loading  head,  than  those  obtained  using  the 
restrictive  head. 

The  above  experimental  trends  did  not  tally  with 
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the  conclusion  of  Gjelsvik  and  Bodner  that  "...for  large  X 
buckling  is  governed  by  the  anti- symmetri cal  modes."  The 
shape  acouired  by  the  arch  is  irrefutably  asymmetrical,  yet 
they  state  this  phenomenon  could  be  observed  only  by  care¬ 
ful  experimentation  as  in  their  own  case  they  witnessed 
very  little.  This  was  due  to  the  very  shallow  arches  they 
investigated.  The  author  assumes  therefore  that  the  above 
phrase  refers  to  the  load  values  and  the  shape  of  the 
resulting  curve. 

Presumably  they  drew  such  a  conclusion  from 
Figure  26,  reproduced  from  their  paper.  This  writer  is  at 
a  loss  to  interpret  their  symmetrical  mode  curve  for  it 
can  infer  only  that  at  the  critical  value,  Py  and  P^  begin 

to  diverge  rapidly.  Gjelsvik  and  Bodner  presented  only  one 
full  set  of  curves,  for  X  =  11.62,  and  this  is  also  reprod¬ 
uced  in  Figure  27;  this  curve,  it  is  claimed,  supports  the 
view  that  after  the  critical  value  of  X  =  10.6,  the  transi¬ 
tional  mode  occurs  before  the  upper  buckling  load  is  reached. 
Hence  it  follows  that  a«p  is  less  than  ay,  and  supposedly 
for  large  values  of  X,  the  situation  given  in  Figure  2$ 
would  arise. 

Such  a  curve  was  not  found  for  any  of  the  experi¬ 
mental  graphs  obtained  by  the  author.  Further,  for  no  value 
of  X  did  the  author,  in  applying  the  G jelsvik-Bodner  theory 
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find  that  Sy  was  greater  than  a^.  Although  a  passing 

glance  at  Figure  29  would  perhaps  indicate  that  P<j> 

* 

occurred  before  Py ,  Figure  30  leaves  no  doubt  that  the 
upper  buckling  load  is  achieved  first  at  deflection  ay. 
Thence  at  deflection  ax  the  secondary  buckling  mode  sets 
in. 

The  critical  value  of  10.6  is  not  derived 
mathematically  in  Gjelsvik  and  BodnerTs  paper,  and  its 
exactitude  is  puzzling.  Taking  the  two  term  series  for 
the  shape  function  (27),  expression  (4 6)  for  the  transi¬ 
tional  buckling  load  is 


P* 

r  t  ss  7T 

12X 


1  + 


0.2245  -  0. !&47rA 
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Manipulation  of  the  portion  under  the  radical 
reveals  that  the  lowest  value  of  X  for  which  an  arch  can 
experience  a  transitional  mode  is  X  =  3.96.  Since  at 

this  value  *  ^  2 

FT  *  PE  =  -IL 
12X 

it  may  be  deduced  that  P-p  occurs  first  at  the  intermediate 
buckling  load,  with  the  mode  reverting  to  symmetry  immed¬ 
iately  upon  further  deflection.  For  subseouent  increases 
in  X,  the  relative  position  of  secondary  buckling  occur¬ 
rence  rises  until  the  mode  occurs  consistently  at  a  =0.2211, 
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and  vanishes  at  a  constant  deflection  value,  just  before 
the  lower  buckling  load  is  reached. 

Employing  a  shape  function  with  six  terms, 
Gjelsvik  and  Bodner  derived  that 

l 

0.16S5  -  CU843n41 " 1 

X2  J  J 

Since  the  lowest  value  of  X  for  which  Pp  is  real  is  now 
X  =  10.3,  and  at  this  value,  secondary  buckling  deform¬ 
ations  are  instantaneous,  it  would  seem  most  unreasonable 
to  surmise  that  at  X  =  10.6,  the  transitional  mode  becomes 
dominant.  The  author,  therefore,  on  the  basis  of  Gjelsvik 
and  Bodner1 s  own  theory,  does  not  agree  with  their  conclu¬ 
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BODNER'S  T  >£OR'.  ,  FOP  PiECES  SI-S4 
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FIG.  I  3  THEORETICAL  AND  EXPERIMENTAL  LOAD  DEFLECTION 
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FIG.  15  T;  EORETICAL  AND  EXPERIMENTAL  LOAD  DEFLECTION 
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FIG.  16  THEORETICAL  AND  EXPERIMENTAL  LOAD  DEFLECTION 
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FIG„  17  THEORETICAL  UPPER  BUCKLING  LOADS 
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FIG.  13  EXPERIMENTAL  LOAD  DEFLECTION  CURVES 
COMPARED  WITH  WIJNGAARDEN  THEORY, 

^  CONSTANT 
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FIG.  19  EXPERIMENTAL  LOAD  DEFLECTION  CURVES 
COMPARED  WITH  WIJNGAARDEN  THEORY, 
y  CONSTANT 
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FIG.  21  EXPERIMENTAL  LOAD  DEFLECTION  CURVES, 
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PIECE  S6 
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COMPARED  TO  CURVE  OBTAINED  3Y 


VAN  WIJNGAARDEN  THEORY 
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FIG.  23  7  'EORE1  ICAL  BUCKLING  LOADS  EXPRESSED  AS 


FIG.  24  PIECL  S7  BUCKLED  UNDER  OWN  WEIGHT 
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FIG.  25  THEORETICAL  AND  EXPERIMENTAL  LOAD  DEFLECTION 
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MG,  26  COMPARISON  WITH  GJELSVIK  AND  BODNER'S  THEORETICAL  BUCKLING  LOAD  CALCULATIONS 
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Flo.  27  COMPARISON  WITH  GJELSVIK  AND  BODNER'S  RESULTS 
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FIG.  28  LOAD  DEFLECTION  CURVE  PREDICTED  BY  GJELSVIK 


AND  BODNER  FOR  LARGE  X  VALUES 


71 


FIG.  29  THEORETICAL  UPPER  AND  TRANSITIONAL  BUCKLING  LOADS 
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CHAPTER  5 
CONCLUSIONS 

For  the  common  consideration  of  symmetrical 
buckling  the  G jelsvik-Bodner  theory  and  the  van  Wijngaarden 
theory  agree  closely  with  each  other,  provided  X  is  greater 
than  20.  It  would  thus  appear  that  the  two  term  series 
solution  is  indeed  most  adequate  to  describe  the  buckled 
shape. 

It  is  difficult  to  assess  why  the  variation  in 
dimensional  elimination  techniques  exists  and  why 
van  Wijngaarden  seems  to  have  used  a  stronger  set  of 
factors  than  Gjelsvik  and  Bodner,  although  both  are 
correct.  Buncai}  in  his  text  on  the  subject^,  indicates 
that  commonly  occurring  quantities,  such  as  flexural 
rigidity,  should  be  used  in  preference  to  an  arbitrary 
selection  of  indices.  The  author1 s  findings  confirm  this. 
Analyses  which  do  not  make  any  such  consideration  initially 
terminate  with  several  choices  for  the  various  indices, 
as  the  existing  boundary  conditions  do  not  supply  enough 
equations  to  satisfy  the  indice  unknowns.  It  appears  that 
in  the  G jelsvik-Bodner  theory,  the  factors  chosen  were 
selected  rather  for  their  arithmetical  convenience,  than 
for  their  relative  non-dimensional  advantage. 
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Both  theories  neglect  any  weight  consideration. 

It  is  felt  that  this  omission  is  a  large  source  of  error 

for  the  arch-springs  considered  in  this  thesis,  and  in  fact 

for  any  very  light  arches  with  a  large  radius  of  curvature. 

In  the  range  Gjelsvik  and  Bodner  investigated,  a  weightless 

consideration  would  be  satisfactory  because  of  the  large 

buckling  load  to  weight  ratios.  However,  as  it  stands,  their 

theory  would  indicate  that  it  may  be  applied  to  any  value  of 

* 

X,  but  since  they  appear  to  suggest  that  Py  attains  large 
proportions  after  the  critical  value  of  10.6,  then  it  would 
follow,  perhaps,  that  such  a  restriction  would  be  overworked. 
Still,  as  far  as  can  be  judged,  their  curve  for  transitional 
buckling  concurs  with  the  author* s  and  if  so,  such  a  weight 
neglection  would  bring  in  errors. 

Theoretical  predictions  were  therefore  higher 
than  results  obtained  by  experimental  investigation,  although 
the  symmetrical  buckling  deflections  recorded  were  in  close 
agreement.  It  would  seemingly  be  beneficial  if  these 
theories  did  include  terms  for  distributed  weight,  albeit 
that  this  would  complicate  considerably  the  succeeding 
derivation  of  the  load  deflection  relationship. 

The  author  disagrees  partially  with  Gjelsvik  and 
Bodner  in  their  treatment  of  the  nature  of  the  asymmetrical 
mode.  He  would  suggest  instead,  that  under  completely 
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static  conditions,  an  arch  will  continue  to  buckle  at  the 
upper  buckling  load ,irregardless  of  its  X  ratio.  Then  if 
X  is  of  sufficient  magnitude,  shortly  after  this  happening, 
the  arch  will  deform  asymmetrically.  With  the  approach  of  the 
second  phase  of  stability  the  asymmetry  will  decrease,  and 
symmetry  will  be  assumed  just  before  the  lower  buckling 
load  is  attained.  Their  postulates  that  transitional 
buckling  distinguishes  itself  by  a  straight  line  plot  on  a 
load  deflection  graph,  and,  that  this  line  intersects  with 
the  symmetrical  curve  at  the  intermediate  load  value,  are  in 
keeping  with  the  findings  of  this  writer.  Further  on  the 
evidence  of  Figures  23  and  26,  he  supports  the  view  that  for 
values  of  X  less  than  approximately  3.5,  buckling  does  not 
occur. 

As  Figure  23  indicates,  the  van  Wijngaarden  theory 
should  be  applied  with  caution  for  low  X  values.  Since  this 
theory  was  derived  specifically  for  T,thin  rings’1,  then  such 
a  limitation  is  to  be  expected,  but  as  a  consequence  the 
theory  gives  no  indication  of  the  initially  changing  nature 
of  buckling. 
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APPENDIX  A 

Young’s  Modulus  for  the  banding  steel  used  in 
the  experimental  work  was  determined  by  employing  a 
test-piece,  which  was  about  a  foot  in  length.  Along 
lightly  scribed  centre  lines,  two  Budd  Metafilm  strain 
gauges  were  affixed  on  opposite  sides  of  the  banding  steel, 
and  a  similar  strain  gauge  was  also  cemented  to  a  dummy 
strip  of  the  banding  steel.  Before  adhesion  of  each 
gauge,  the  surface  of  the  steel  was  rubbed  with  fine 
sandpaper  and  cleaned  with  trichloroethylene. 

The  test  length  was  then  placed  in  a  tension 
tester,  and  it  was  wired  to  a  strain  indicator  and  a 
switching  box.  Load  was  applied  on  the  strip  several 
times  before  readings  were  noted  and  the  results 
obtained  are  shown  graphically  in  Figure  31. 

Young’s  Modulus  was  then  calculated  as  follows: 
Area  of  specimen  *=  0.75  x  .0150 

»  0.1125  sq.  inches 

E  =  P  1 
e  0.1125 

Taking  the  value  of  P  from  the  gradient  of  the  line  in 

e 

Figure  31,  then 

E  »  269  1 

800  0.1125 

E  =  29.9  x  10^  p.s.i. 


so  that 


. 

,  • 

, 

* 

- 

. 


300 


78 


F3G„  31  LOAD  STRAIN  CURVE  TO  DETERMINE  YOUNG’S  MODULUS 


APPENDIX  B 


IBM  1620  Computer  programme  yielding 
theoretical  load  deflection  curves  for  varying 
X  values. 
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